Weakly Pinned Bose Glass vs. Mott Insulator Phase in Superconductors 
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We study the properties of the Bose glass phase of localized flux lines in irradiated superconductors 
near the matching field _B$. Repulsive vortex interactions destroy the Mott insulator phase predicted 
to occur at B — B§. For ratios of the penetration depth to average defect distance X/d < 1 remnants 
of the Mott insulator singularities remain visible in the magnetization, the bulk modulus, and the 
magnetization relaxation, as B is varied near B$. For A > d, the ensuing weakly pinned Bose glass 
is characterized by a soft Coulomb gap in the distribution of pinning energies. 

PACS numbers: 74.60. Ge, 05.60,+w 



For the application of high-T c (type II) superconduc- 
tors in external magnetic fields, an effective flux pin- 
ning mechanism is essential in order to minimize the re- 
sistive losses through Lorentz-force induced vortex mo- 
tion. Specifically, columnar defects, i.e., linear damage 
tracks in the material caused by heavy-ion irradiation, 
have emerged as very effective pinning centers For 
such systems, a continuous vortex localization transi- 
tion at Tbg from an entangled flux liquid to a disorder- 
dominated Bose glass phase was predicted [§,0 , and sub- 
sequently found in experiment jO. In addition, it was 
suggested that when the vortex and defect densities are 
equal, each flux line would be attached to one pin, leading 
to a Mott insulator phase within the Bose glass j| . Re- 
cent measurements of the magnetization relaxation rate 
at low temperatures ||,|5| , and of the reversible magneti- 
zation itself 1^,0, have been interpreted as signatures of 
this Mott insulator, and inspired further theoretical in- 
vestigations ^-H- Yet so far, the influence of the repul- 
sive vortex interactions, whose range is set by the London 
penetration depth A, has not been carefully studied. 

The aim of this Letter is to explore numerically how 
the repulsive vortex forces affect the predicted Mott insu- 
lator as the magnetic field B is varied near the matching 
field Using a random defect distribution with aver- 
age defect distance d, we extend earlier work which was 
limited to B < 5$ Q. Our findings are: 1) At _B$ the 
Mott insulator, which is characterized by a hard gap A in 
the distribution of pinning energies 17(e) near the chem- 
ical potential /i, exists only for extremely short-range 
interactions X/d — > (Fig. |l|, left inset); with increasing 
X < d the gap quickly fills and, therefore, the Mott insu- 
lator is destroyed (Fig. [|, right top inset). For long-range 
interactions A > d a soft Coulomb gap Q described by 
g(e) oc |e — ^| Scff emerges with an effective gap expo- 
nent s ff (Fig. [j], bottom inset). 2) By measuring the 
reversible magnetization and the effective IV exponent, 
which is related to the magnetization relaxation rate, we 
are able to explain experimental results as remnants 
of the Mott insulator for A < d, but not as true signatures 
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FIG. 1. SBG / WBG crossover line (see text) vs. X/d. The 
insets show the pinning energy distribution g(e) vs. e — fj, for 
X/d -> (left), X/d = 1/4 (right top) and X/d = 5, all at B<s>. 



of this distinct thermodynamic phase itself. 3) We can 
identify a crossover line (Fig. [j]), discriminating between 
the strong Bose glass (SBG, with all vortices localized 
by defects) and the weak Bose glass (WBG, with vortex 
bundle pinning) ||,|^]. 

The theoretical description of the Bose glass is based 
on the following free energy of Ny flux lines, described 
by their 2d trajectories Yi(z) as they traverse a sample of 
thickness L in a mag netic field B || z [§,§, 
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This consists first of an elastic line tension term, with tilt 
modulus e\ . The second term denotes the interaction en- 
ergy of all vortex pairs (local in z), where nj = \r{ — Tj\ 
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and V(r) = 2eoKo(r / X) is the screened repulsive vor- 
tex potential, with the modified Bessel function Kq(x) oc 
— log(x) as x — > 0, and Kq(x) oc a; -1 / 2 exp(— x) for x — > 
oo. The energy scale is set by eo = (tfio/AnX) 2 . The last 
term describes Nd columnar pins (|| B), modeled by z- 
independent square well potentials Vb with average spac- 
ing d, centered on randomly distributed positions {R/c}. 
The typical defect radius is Ck ~ 50A, and related to the 
pinning strengths U k via U k w (e /2) log[l + (c k /V2£) 2 } 
The ion beam dispersion induces a distribution of 
the Uk with width w = \J (5U 2 ). 

The mathematical analysis of Eq. (|l|) exploits a map- 
ping of the statistical mechanics of this free energy 
of directed lines to the quantum mechanics of two- 
dimensional bosons subject to point disorder ||||. Here, 
we focus on the ground states properties of Eq. (Q) for 
B ~ B<$ at zero temperature, which is a fair approxima- 
tion for the Bose glass phase for T < T\ « 0.6 . . . 0.8 Teg- 
In this regime, thermal wandering of vortices can be ne- 
glected (|]ll],[l3|,|l4| i an d the vortices will essentially be 
straight which allows us to ignore the elastic energy. This 
leaves us with a static two-dimensional problem of Ny 
interacting "particles" and Nd defects. For computa- 
tional reasons we represent our problem on an underly- 
ing triangular grid with N "lattice" sites and therefore 
N — Nd interdefect sites ( "interstitials" ) . The effective 
Hamiltonian then reads 

, N No 

HeB = o X] n i n jV( r ij) + X! Uktk ' ( 2 ) 
i?j fc=l 

where {rii — 0, 1} represent the site occupation numbers. 
If we take £ = 10A and cq = 50A, we obtain for the 
"bare" defect pinning energies t k = —{Uk) + Wf., with 
(Uk) — 0.65 and width w = 0.1 in units of 2eo. For sim- 
plicity we assume a flat distribution of pinning energies 
around its mean, P(wk) = Q(w — \wk\)/2w, where 6 de- 
notes the step function. The simulations are performed 
by randomly distributing Ny vortices and Nd defects 
on the triangular grid and minimizing the total energy 
with respect to single-particle transfers, thereby obtain- 
ing pseudo-groundstates for the Hamiltonian (||) f p"5"[ . 
Our simulations are mostly carried out with N = 1600, 
Nd = N/16 =100, using periodic boundary conditions. 
We study physical quantities as functions of the inter- 
action range X/d = 1/4, 1/2, 1, 5, and the filling fraction 
/ = Ny /Nd = B/B&, in the interval 0.2 < f < 3. All 
distances are given in units of the lattice constant a of 
the triangular grid. The lower limit for the interaction 
range is determined by A = a, since even smaller values 
effectively correspond to A — ► 0. We typically take an 
average over 50 different realizations of the disorder. 

We have compared our findings for / < 0.6 with those 
obtained in Ref. where it was assumed that all vor- 
tices remain bound to defects, irrespective of the value of 
A, and therefore continuously spaced random positions 
could be used. Indeed, we were able to reproduce the 



previous results for filling fractions / = 0.2 and / = 0.4 
quantitatively. We have also tried to render the artifi- 
cial grid finer by keeping Nd constant and increase N 
from 1600 to 3600, the largest sample we could study, 
and did not detect any lattice dependence for both grid 
sizes. Hence, we believe that our results should be largely 
insensitive to the underlying lattice representation, and 
provide a fair approximation to a more realistic contin- 
uum description. Also, we could not find any significant 
finite-size effects when keeping N — 3600 and X/d fixed 
and increasing Nd from 36 to 100 and 225. 

A natural first question to ask is how many flux lines 
are depinned as a result of the vortex interactions, de- 
pending on the values of / and X/d. Within the Bose 
glass one may discriminate between the SBG and the 
WBG, where the latter is characterized by a markedly 
reduced localization temperature Tbg and critical cur- 
rent J c ||]. For X/d -C 1, the crossover between these 
regimes is expected to occur for B ~ _B$ (§]. Once inter- 
actions become strong, the WBG will appear well below 
B<p. Fig. |l| shows the filling fraction / occ at which 10 % 
of the vortices are depinned, as a function of X/d, which 
we tentatively take as a criterion to define a crossover 
line between the SBG and WBG. One observes a rather 
strong dependence of this line on X/d, which is shifted 
well below B$ as soon as A w d. Only for X/d < 1/4 
does the line remain above B§; yet obviously this is a 
prerequisite for the Mott insulator phase to appear. 

In order to see if the Mott insulator persists for a finite 
interaction range and a random distribution of pinning 
sites, we compute the single-particle density of states 
g(e), i.e. the distribution of (interacting !) pinning en- 
ergies, where = J2j^i n jV( r ij) + F° r long-range 
interactions (A 3> d) one expects a soft "Coulomb" gap 
in g[e) at the chemical potential e = fx, separating oc- 
cupied (e < fi) from empty (e > fi) states; near u, g(e) 
should vanish according to g(e) oc |e — ^| Scff |jl^,0. In 
the limit X/d — > 0, on the other hand, the Mott insula- 
tor phase at / — 1 is characterized by the appearance 
of a hard gap separating the occupied defect states at 
e = —(Uk) ± w from the unoccupied states at e = 
(Fig. |l|, left inset). The right insets in Fig. (|l|) show g(e) 
at / = 1 for X/d = 1/4 (top) and X/d — 5 (bottom). 
The bottom figure exhibits a wide Coulomb gap, which 
we also find (though narrower) for X/d = 1. 

The top inset, however, only shows a fairly flat den- 
sity of states for the occupied sites, which is somewhat 
depleted for e ~ /i, and rises sharply for e > fi. The fact 
that all states below /i are continuously filled implies that 
even rather short-range interactions suffice to destabilize 
a distinct thermodynamic Mott insulator phase ||. We 
have also looked at X/d = 1/8 by increasing N to 3600, 
using only Nd — 56, whence d — \JN/Nd ~ 8, and 
even here we still found some states in the vicinity of fi. 
We are therefore led to conclude that for the true Mott 
insulator to appear in this system with a random spa- 
tial distribution of pinning sites, a necessary condition is 
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FIG. 2. Log-linear plots of the bulk modulus en (top) and 
the magnetization M (bottom) vs. / for \/d = 1/4 (left scale 
in the bottom plot) and X/d = 1/2 (right scale). 



A -C d. We speculate that for a more regular array of 
columnar defects fl6|] , the Mott insulator at T = may 
persist to considerably larger interaction ranges || . 

Next, we have determined the chemical potential /j, 
as a function of /, i.e. the H(B) curve, defining fi = 
(e mn + e m ; n )/2 with e max the maximum site energy of 
the occupied states, and e m ; n the minimum of the un- 
occupied site energies, and then averaging over disorder 
]l2fl . From this we obtain the bulk modulus of the sys- 
tem via en = dfi/df = Nod^i/dNv- A divergence of c u 
at B$ would indicate the appearance of the Mott insula- 
tor, which is incompressible (adding another vortex costs 
an energy A) and in this respect similar to the Meissner 
state. Correspondingly, its signature should be a sharp 
jump in fi at / = 1 [for A — > the height of this jump 
is given by the hard gap A = (Uk) in g(e)]. Fig. || (top) 
shows that while the divergence is smeared out already 
for short-range interactions (X/d = 1/4), yet a quite 
pronounced peak in en remains, shifted downwards to 
/ » 0.85. For larger interaction range (X/d = 1/2) the 
jump of fi is displaced to even smaller /(ss 0.55) and 
less marked; for A > d, the bulk modulus becomes essen- 
tially constant over the entire range of /. In both cases 
depicted here, the location of the peak in cu coincides 
with that value of / at which a sizeable number of vor- 
tices actually leave the defect sites. The fact that the 
average bulk modulus is enhanced for larger A is simply 
an energetic effect, since a vortex entering the system has 
to overcome a higher energy barrier due to the interac- 
tions. The behavior of cu confirms again that no true 
Mott insulator exists even for relatively short-range in- 
teractions X/d — 1/4 or 1/2. We do observe, though, a 
distinctive "lock-in" structure for X/d < 1, which com- 
pletely disappears for X/d > 1. 

In order to make further contact with experiments, 



we measured the total energy G, which yields the re- 
versible magnetization via the thermodynamic relation 
M = -dG/dB oc -Ny X dG/df. Fig. | (bottom) depicts 
M as a function of / on a log-linear plot for X/d = 1/4 
and X/d = 1/2. The data for X/d = 1/4 show a pro- 
nounced minimum in M at / rs 1, embedded in a slow 
logarithmic growth as / increases. The second plot for 
X/d = 1/2 hardly displays any structure aside from a 
shallow dip near / « 0.6. This feature is completely 
absent for X/d > 1, where only the log(v / /) increase re- 
mains, resembling the magnetization curve of an unirra- 
diated superconductor PJt|JTo||. The observed behavior 
of M at X/d — 1/4 qualitatively agrees very well with re- 
cent measurements performed by van der Beek et at in 
an irradiated BSCCO crystal §, who find (at T » T x ) 
a pronounced dip in M centered near B = B$. The 
disappearance of this minimum in the experiments as 
T — ► TgQ may be at least partially due to the increase 
of the London penetration depth X(T), in addition to 
the entropic renormalizations studied in Ref. [fTo|| . Our 
simulations cannot explain, however, the magnetization 
minima found at / > 1 in BSCCO tapes |Q. 

We would like to point out that the maximum in en 
occurs at lower values of / than the minimum in M 
(compare, e.g., our data for X/d — 1/4), which can 
be understood as follows. Starting from the relation 
M cx B — H, the bulk modulus may be written as 
c u cx dH/dB cx const. - dM/dB = const. + d 2 G/dB 2 ; 
thus the maximum of en occurs at the location of the 
steepest negative slope in M(B), which has to be at a 
smaller B than the minimum in M. Only in the "true" 
Mott insulator phase will the singularities both in en and 
M coincide at B = B$. 

We now briefly return to the case of long-range inter- 
actions X/d > 1, where a soft Coulomb gap appears in the 
distribution of pinning energies (r.h.s. inset in Fig. 0). 
Notice that the energetically favorable sites include both 
the defect and high-symmetry interstitial positions; thus, 
the original spatial randomness is in effect smoothencd, 
and the gap exponent s e g increased as compared to a sit- 
uation where no interstitials are allowed At large 
/, a substantial fraction of the vortices will leave the 
pin positions, and as a consequence, Bose glass behavior 
can only be seen at low temperatures Yet, then the 
pinning via the repulsive forces exerted by neighbors that 
are still attached to columnar defects is very effective; we 
find that the potential minima are roughly equally deep 
and wide for all vortices in a given pseudo-ground state 
configuration, irrespective of the flux line occupying a de- 
fect or an interstitial site. This gives us confidence that 
the single-vortex density of states may indeed be used 
to infer low-current transport properties in the variable- 
range hopping regime 

In the localized Bose glass phase, vortex transport at 
low currents J <C J c is expected to occur via variable- 
range hopping [H, in analogy with doped semiconduc- 
tors Il2| . In the spirit of the thermally-assisted flux flow 
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FIG. 3. Transport exponent l/p e g vs. / for X/d = 1/2, 1, 5. 



model, this leads to a highly nonlinear IV characteristics 
£ = po J exp[— Db(J)Ab7]i w hh effective (free) energy 
barriers that diverge according to Ub(J) = Uo (Jo/ ' J) p " s 
as J — > 0. For short-range interactions, p c $ should be 
given by the 2d Mott variable-range hopping exponent 
Po = 1/3 ||. Yet, once the vortex interactions become 
long-range, the emergence of a Coulomb gap in the dis- 
tribution of pinning energies near [i leads to a consid- 
erable enhancement of flux pinning with effective expo- 
nents up to p c ff w 0.7, if / = 0.1 JyJ. Such values for 
p e ff in the range between 1/3 and 1 (as suggestive of the 
vortex half-loop excitations dominating for intermediate 
currents PJTl[|) were found in recent magnetization re- 
laxation experiments extending up to -B$ jH]. 

We have calculated p c g from the IV characteristics ob- 
tained by integrating over g(e) for e > \x (see Ref. @). 
In Fig. |3] we plot p~g as a function of / for various val- 
ues of X/d. Remarkably, p e f{ increases with growing vor- 
tex density for / < 1, as opposed to the earlier simu- 
lations where no interstitial positions were allowed and 
to the contrary a decrease with / was found, owing to 
the fact that the system had to accommodate with the 
underlying randomness [j]~T]| . Now, however, with favor- 
able interstitial sites being available, the disorder effects 
are screened by the interactions as the vortex density in- 
creases. This leads to stronger correlations, and in fact 
more mean-field like behavior with larger values of p c s 
as is apparently seen in experiment Q]. For X/d = 1/2, 
we find a marked minimum in p~^ at/« 1.4 (see Fig. |) 
obviously due to a delicate interplay of disorder and 
correlation effects. Assuming that relaxation times are 
determined by the variable-range hopping energy bar- 
riers (vortex superkinks), log(t/to) ~ UB{J)/kBT |Q, 
one finds for the magnetization relaxation rate S = 
-dM/dlogt = -dlogJ/dlogt = k B T/p cB U{J); there- 
fore S(B) cx l/p e s(f). Hence, the minima in S(B) 
near B « 1.4i?$ detected by low-temperature relaxation 



experiments on YBCO B, and on a Tl-compound [BJ, 
might be explained by the maximum in p e g we find with 
our simulations. At elevated temperatures, this structure 
should disappear quickly in this WBG regime. 

In summary, we have investigated the Bose glass phase 
of flux lines localized by columnar defects near the match- 
ing field, taking the vortex interactions properly into ac- 
count. For a random spatial distribution of pinning cen- 
ters, the proposed Mott insulator phase is destabilized 
by the repulsive forces between the flux lines, but for a 
moderate interaction range X/d < 1 interesting "lock-in" 
effects remain visible in the reversible magnetization, the 
bulk modulus, and the magnetization relaxation, which 
at least qualitatively explains a number of recent exper- 
iments. For larger X/d > 1, the distribution of pinning 
energies displays a wide Coulomb gap, and sites available 
for variable-range hopping processes include both defect 
as well as high-symmetry interstitial positions. 
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